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Abstract 

The latest results of the ATLAS and CMS experiments point to a preferred narrow Higgs mass 
range (m^ ~ 124 — 126 GeV) in which the effective potential of the Standard Model (SM) develops 
a vacuum instability at a scale 10 9 — 10 11 GeV, with the precise scale depending on the precise 
value of the top quark mass and the strong coupling constant. Motivated by this experimental 
situation, we present here a detailed investigation about the stability of the SM ++ vacuum, which 
is characterized by a simple extension of the SM obtained by adding to the scalar sector a complex 
SU(2) singlet that has the quantum numbers of the right-handed neutrino, H", and to the gauge 
sector an U(l) that is broken by the vacuum expectation value of H" . We derive the complete set 
of renormalization group equations at one loop. We then pursue a numerical study of the system to 
determine the triviality and vacuum stability bounds, using a scan of 10 8 random set of points to 
fix the initial conditions. We show that, if there is no mixing in the scalar sector, the top Yukawa 
coupling drives the quartic Higgs coupling to negative values in the ultraviolet and, as for the SM, 
the effective potential develops an instability below the Planck scale. However, for a mixing angle 
0.05 < a < 0.35, with the new scalar mass in the range 500 GeV < mi» < 5 TeV, the SM ++ 
ground state can be absolutely stable up to the Planck scale. These results are largely independent 
of TeV-scale free parameters in the model: the mass of the non-anomalous U(l) gauge boson and 
its branching fractions. 



*On leave of absence from CPHT Ecole Polytechnique, F-91128, Palaiseau Cedex. 



I. INTRODUCTION 



The CERN Large Hadron Collider (LHC) has begun a bound and determined exploration 
of the electroweak scale. Recently, the ATLAS [1] and CMS [2] collaborations presented an 
update of the Higgs searches, independently combining about 5 fb _1 of data collected at 
y/s = 7 TeV and more than 5 fb _1 at y/s = 8 TeV. The excess at 125 GeV that was evident 
already in data from the 7 TeV run [3, 4] has been consistently observed by both experiments 
in the 77 invariant mass spectrum with a local significance of 4.5a and 4.1a, respectively. 
In addition, an excess of 4 leptons events (with rriu ~ 125 GeV) which can be interpreted 
as a signal of the h — > ZZ* — > 4£ decay, is observed by both experiments with a significance 
of 3.4a and 3.2a, respectively. The CMS experiment also presented updated Higss boson 
searches in W + W~ (a broad excess in the invariant mass distribution of 1.5a is observed), bb 
(no excess is observed), and rf (no excess is observed) channels. More recently, the ATLAS 
Collaboration reported a 2.8a deviation in the h — > W + W~ — > 2£u decay channel [5]. When 
combining the data from the 7 TeV and 8 TeV runs, both experiments separately have 
reached the sensitivity to the new boson with a local significance of 5a [6, 7]. Very recently 
recently, the CDF and DO collaborations published an update on searches for the Higgs 
boson decaying into bb pairs, using 9.7 fb _1 of data collected at y/s = 1.96 TeV [8]. They 
reported a 3.3a deviation with respect to the background-only hypothesis in the mass range 
between 120 - 135 GeV. 

All in all, LHC data strongly suggest that the observed state feed the electroweak sym- 
metry breaking, and is likely the Higgs boson. However, it remains to be seen whether its 
trademarks, particularly the production cross section and decay branching fractions, agree 
with the very precise prediction of the Standard Model (SM). The decay channels that are 
most sensitive to new physics are the loop induced decays h — > 77 and h — > jZ. Interest- 
ingly, the most recent analyses [9, 10] of the combined LHC data seem to indicate that there 
is a deviation from SM expectations in the diphoton channel at the 2.0a — 2.3a level; see, 
however, [11]. 

From a theoretical perspective some modification of the Higgs sector has long been ex- 
pected, since the major motivation for physics beyond the SM is aimed at resolving the 
hierarchy problem. Even if one abandons such a motivation for new physics there are still 
enduring concerns about the stability of the electroweak vacuum, which have been exacer- 
bated by the new LHC data that points to rrih — 125 GeV. 

Next-to-leading order (NLO) constraints on SM vacuum stability based on two-loop renor- 
malization group (RG) equations, one-loop threshold corrections at the electroweak scale 
(possibly improved with two- loop terms in the case of pure QCD corrections), and one-loop 
effective potential seem to indicate rrih ~ 125 — 126 GeV staurates the minimum value that 
ensures a vanishing Higgs quartic coupling around the Planck scale (Mpi), see e.g. [12- 
22, 39]. However, the devil is in the details, a more recent NNLO analysis [23] yields a very 
restrictive condition of absolute stability up to the Planck scale 



m h > 



GeV. (1) 



When combining in quadrature the theoretical uncertainty with experimental errors on the 
mass of the top (mt) and the strong coupling constant (a s ), one obtains rrih > 129±1.8 GeV. 
The vacuum stability of the SM up to the Planck scale is excluded at 2a (98% C.L. one 
sided) for rrih < 126 GeV [23]. Achieving the stability will neccesarily impose constrains on 
physics beyond the SM. 
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Very recently we have put forward a (string based) Standard-like Model [24]. Motivated 
by the above, here we study the vacuum stability of its scalar sector. The layout of the 
paper is as follows. In Sec. II we briefly review the generalities of our model and we derive 
the RG equations. In Sec. Ill we present our results and conclusions. 

II. RG EVOLUTION EQUATIONS OF SM++ 

Very recently, we engineered the minimal extension of the SM that can be embedded into 
a Superstring Theory endowed with a high mass string scale, M s < Mpi [24]. The gauge 
extended sector, f/(3)sX SU (2) L x U(l)j R x U has two additional U (1) symmetries and 
thus we refer to our model as SM ++ . The origin of this model is founded on the D-brane 
structure of string compactifications, with all six extra dimensions O(M^) [25-27]. The low 
energy remnants of the D-brane structure are the gauge bosons and Weyl fermions living 
at the brane intersections of a particular 4-stack quiver configuration [28]. A schematic 
representation of the D-brane construct is shown in Fig. 1. The general properties of the 
chiral spectrum are summarized in Table I. 

The resulting U(l) content gauges the baryon number B [with U(1)b C U(3)b], the 
lepton number L, and a third additional abelian charge Ir which acts as the third isospin 
component of an SU(2) R . Contact with gauge structures at TeV energies is achieved by a 
field rotation to couple diagonally to hypercharge V M . Two of the Euler angles are determined 
by this rotation and the third one is chosen so that one of the U(l) gauge bosons couples 
only to an anomaly free linear combination of Ir and B — L. Of the three original abelian 
couplings, the baryon number coupling g' 3 is fixed to be \/l/Q of the QCD coupling <? 3 at the 
string scale. The orthogonal nature of the rotation imposes one additional constraint on the 
remaining couplings g[ and g' 4 [29]. Since one of the two extra gauge bosons is coupled to an 
anomalous current, its mass is 0(M S ), as generated through some Stiickelberg mechanism. 1 
The other gauge boson is coupled to an anomaly free current and therefore (under certain 
topological conditions) it can remain massless and grow a TeV-scale mass through ordinary 
Higgs mechanisms [31]. 

Electroweak symmetry breaking is achieved through the standard Higgs doublet H. The 
spontaneous symmetry breaking of the extra non-anomalous U(l) is attained through an 
577(2) singlet scalar field H", which carries L and Ir numbers, and acquires a vacuum ex- 
pectation value (VEV) at the TeV scale. With the charge assignments of Table I there are 
no dimension 4 operators involving H" that contribute to the Yukawa Lagrangian. This 
is very important since H" carries the quantum numbers of right-handed neutrino and its 
VEV breaks lepton number. However, this breaking can affect only higher-dimensional op- 
erators which are suppressed by the high string scale, and thus there is no phenomeno logical 
problem with experimental constraints for M s higher than ~ 10 14 GeV. Herein we remain 
agnostic with respect to supersymmetry (SUSY) breaking and the the details of the low 
energy effective potential. However, we do subject the choice of quantun numbers for H" to 
the stringent holonomic constraints of the superpotential at the string scale. This forbids 



A point worth noting at this juncture: SM can also be embedded in a 3-stack quiver comprising (only) 
one additional U(l) symmetry, U(3) x SU(2) x U(l) [30]. The extra gauge boson is anomalous and must 
grow a Stiickelberg mass <~ M s . In this D-brane model the running of the quartic Higgs coupling would 
reveal a vacuum instability around 10 10 GeV [23]. 
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FIG. 1: Pictorial representation of the U(3)b x SU(2)l x C/(1)l x £7(1)/^ D-brane model. 

TABLE I: Chiral spectrum of SM++. 
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the simultaneous presence of scalar fields and their complex conjugate. As an illustration, 
if the quantum numbers of H" are those of N R , then higher dimensional operators such 
as N R N R H", which can potentially generate a Majorana mass, are absent. Because of 
holonomy this absence cannot be circunvented by including N rN r H"* 2 . 
The scalar Lagrangian of SM ++ is 



& 8 = (2W) T V^H + (V»H") j V^H" - V(H, H" 



(2) 
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where 

V (H, H") = fi 2 \H\ 2 + fi' 2 \H'f + Ax \H\ 4 + A 2 \H"f + A 3 \H\ 2 \H"\ 2 (3) 
is the potential and 

V„ = d„ - ig^Gl - ig' 3 Q B C, - ig 2 r a W; - ig[Q lR B, - ig' 4 Q L X, (4) 

is (in a self-explanatory notation [24]) the covariant derivative in the field basis shown in 
Fig. 1. Next, we impose the positivity conditions [32] 

A!>0, A 2 >0, A 1 A 2 >^A^. (5) 

If the conditions (5) are satisfied, we can proceed to the minimisation of V(H, H") as a 
function of constant VEVs for the two Higgs fields. Making use of gauge invariance, it is 
not restrictive to take 



(H) = ^=( I , A and (H") = ±= (v" + h 2 



») , (6) 



with v and v" real and non-negative. The physically most interesting solutions to the 
minimisation of (3) are obtained for v and v" both non-vanishing 

2 — A 2 /i 2 + |A 3 /i" 2 2 — Ax/i" 2 + |A 3 /i 2 
v = \ \ i A 2 and v =^~ x TT2— • ( 7 ) 

AiA 2 — ^ A 3 AiA 2 — |A 3 

To compute the scalar masses, we must expand the potential (3) around the minima (7). 
We denote by h and h" the scalar fields of definite masses, rrih and rrih" respectively, taking 
m h ~ 125 GeV, v = 246 GeV, and m h „ 3> m^. After a bit of algebra, the explicit expressions 
for the scalar mass eigenvalues and eigenvectors are given by 

ml = \y + X 2V " 2 - (X lV 2 - X 2 v" 2 ) 2 + (X 3 vv") 2 , (8) 

X x v 2 + X 2 v" 2 + ^ (X lV 2 - X 2 v" 2 ) 2 + (X 3 vv") 2 , (9) 

(10) 



m 



h \ _ I cos a — sin a; \ / hi 
h" J \ sin a cos a J \h 2 

where a G [— 7r/2,7r/2] also fulfils 

sin 2a = — = , (11) 

y/(X lV 2 - X 2 v" 2 ) 2 + (A 3 W) 2 

W - \2v" 2 h0 . 
cos2a = — . 12 

^{x lV 2 - x 2 v" 2 y + (x 3 w»y 

Now, it is convenient to invert (8), (9) and (11), to extract the parameters in the Lagrangian 
in terms of the physical quantities rrih, rrih" and sin 2a 
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(13) 
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As we noted above, the fields C^X^B^ are related to Y^Y^Y^' by an Euler rotation 
matrix [33], 

CoCif} —C^S^p + S^SgC^p S<f,S^ + C<pSgC^ \ 

CgS^p C^Cip + S^SgS-ip —S^C^p + C^SgS^p I . (14) 
— Sg S^Cg CcffCg / 

Hence, the covariant derivative for the U(l) fields in Eq. (4) can be rewritten in terms of 
Y^ Y' and Y" as follows 

Dp = df, - iY^ {-Sgg[Q lR + C e S^g' 4 Q L + CeC^g' 3 Q B ) 

~ [CgS ( f > g' 1 Qi R + (C^C^ + SeStpSj/j) g 4 Qi + (C^SgS^ — C^S^g^Qs] (15) 
— iY'p [CoC ( pg' 1 Qi R + (—C^S^ + C^SgS^) g'^Qi + (C^C^Sg + S^S^p) %Qb\ ■ 

Now, by demanding that Y^ has the hypercharge 

Qy = C\Qi R + c 3 Q B + c 4 Q L (16) 

we fix the first column of the rotation matrix R 

YtiCsgy/g's ... \ 

^c 4 ^M ••• h (17) 

YtxCigv/g'i ■■■ J 

and we determine the value of the two associated Euler angles 

= -arcsin[ci5(y/54] (18) 

and 

V> = arcsm[c 4 #y/ (# 4 C<?)] , (19) 

with c\ = 1/2, c 3 = 1/6, C4 = —1/2, B = Qb/3 and L = Ql- The couplings ^ and g' A are 
related through the orthogonality condition, P(gy, g[, #3, #4) = 0, yielding 



C4 \ 1 / C3 \ / Ci 



04/ #r V^3/ \9[ 



(20) 



with g 3 fixed by the relation of U (N) unification 



g' 3 (M s ) = -^=g 3 (M s ). (21) 

Next, by demanding that Y" couples to an anomalous free linear combination of Ir and 
B — L we determine the third Euler angle 

tw * = - s » l e j C s* + fc ■ (22) 

The absence of abelian, mixed, and mixed gauge-gravitational anomalies is ensured utilizing 
the generalized Green-Schwarz mechanism, in which triangle anomalies are cancelled by 
Chern-Simons couplings. In the V-basis we require the (mass) 2 matrix of the anomalous 
sector to be diag(0, M' 2 , 0). For the heavy field we take M' ~ M s and therefore Y^ ~ Z'^ 
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decouples from the low energy physics. The non-anomalous gauge boson, ~ Z" grows a 
TeV-scale mass via H" [24]. 

Altogether, the covariant derivative of the low energy effective theory reads 

V, = d,- ig 3 T a G; - ig 2 r a W; - ig Y Q Y Y, - ig y Q Y Y^ - ig B _ L (B - L) Y[, (23) 

where 

9y — —\Seg'\ — —^CeS^g^ = 7;C 9 C^g' 3 

g y = 2C e C rt> g' 1 (24) 
gB-L = "Sg'ziC^C^Se + S^S^) — CeC ( pg' 1 . 

Finally, a straightforward calculation leads to the RG equations for the five parameters in 
the scalar potential [34] 

,//2 / .,2 



dt 
dfi" 2 



dt 
~dt 



16tt 2 



8A 2 + 4^A 3 -24^_ L 



^ (24A? + XI- 61? + \g\ + jj<£ + ^ + ^ + \gWy + \ 9 y 



+ 12X ± Y 2 - 9\ l9 2 2 - 3X i9y - 3AtfJ) , 

dA 2 1 

"dt" ~ 8^ 

dA 3 A 3 



(25) 



(l0X 2 2 + Xj + 48gjj_ L -24X 2 g 2 B _ L ) , 



^ + 4A 2 + 2A 3 + 3Y t 2 - ^gj - ^g 2 - ^g 2 - \2g 2 B _^j + A ^ ^ 



8tt 2 

where t = In Q and Yj is the top Yukawa coupling, with 



d^ 
dt 



Y /9 



16tt 2 



9 2 _ 17 2 
3 ~~ ^^2 — \2^ Y 



17 2 2 2 



12 



9y 



t^9b-l - 7^9y9B-L 



(26) 



and F/ 0) = 

dg 3 
dt 

dg2 
dt 

dgv 
dt 

dgB-L 
dt 
dgy 
dt 

where n n 



V2m t /v. The RG running of the gauge couplings follow the standard form 



gf 
16tt 2 

jL 

16tt 2 
1 



-11 + 3% 

22 



1M 

16 7T 2 ' 



3 + 3^ + 6 
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16vr 2 

= ^ [A (B - L)( ^ L) ^L L + 2A( fl -^_ Lay + A YY g B _ L g 2 } , 

= ^ [A YY 9 y (g 2 y + 2g 2 Y ) + 2A< b -^b-l(sJ + <£) + ^^^L^] , 

= 3 is the number of generations and 
2 . 1 - 

= A ba = - ^ QajQbJ + o Qa,sQb,s , 



(27) 



(a,b = Y, B-L), 



(28) 



with / and s indicating contribution from fermion and scalar loops, respectively. 
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FIG. 2: An exhibition of the SM ++ vacuum stability patterns in the m/j», a, Mz" parameter 
space. The analysis is based on a scan of 10 8 trial random points with M s = 10 14 GeV. The points 
yielding a stable vacuum solution up to M s are blue-printed, those leading to unstable vacuum 
solutions are red-printed, and points giving runaway solutions are purple-printed. 

III. RESULTS AND CONCLUSIONS 

To ensure perturbativity of g' A between the TeV scale and the string scale we find from 
(20) that g[ > 0.4845. We also take g[ < 1 in order to ensure perturbativity at the string 
scale. Let us first study the region of the parameter space constrained by g[(M s ) ~ 1. The 
string-scale values of the other abelian couplings are fixed by previous considerations (20) 
and (21). The Euler angles at M s are also fixed by (18), (19), and (22). All the couplings 
and angles are therefore determined at all energies through RG running. As an illustration 
we set M s = 10 14 GeV; this leads to g' 3 (M s ) = 0.231, g' 4 (M s ) = 0.232, i/;{M s ) = -1.245, 
0(M 8 ) = -0.217, and <f>(M a ) = -0.0006. Next, we define Q min = 4 TeV and normalize 
t = ln(Q/4 TeV) and t max = ln(A/4 TeV). Finally, we run the couplings and angles down 
to the TeV region: g[ = 0.406, g' 3 = 0.196, g' 4 = 0.218, 9 = -0.466, i/j = -1.215, and 
(f) = -0.0003. 

There are a priori three free parameters to be fixed at the TeV-scale: (v", a, rrih"). The 
initial values of gy, gy and gs-L are then fixed by previous considerations (24). Actually, 
using the relation Mz" = g[C ( j > v" /Co [24], we adopt {Mz», a, nth") as the free parameters 
of the model. 2 For M s = 10 14 GeV, we perform a scan of 10 8 trial random set of points, 
(Mz», a, rrih"), and using (13) we obtain the initial conditions (A]° , \f \ X 3 ) to integrate 
(25). For each set of points, we verify that the positivity condition (5) is fulfiled all the way 
to A = M s . The 10 8 trials are duplicated for M s = 10 16 and M s = 10 19 GeV. Our results 
are encapsulated in Figs. 2 to 9, and along with other aspects of this work are summarized 
in these concluding remarks: 

• In Fig. 2 we show the results for the entire scan at M s = 10 14 GeV. The points yielding 
a stable vacuum solution up to M s are blue-printed, those leading to unstable vacuum 



2 For M s = 10 14 GeV, the v" ^ M z » relation implies that if 7 TeV < v" < 15 TeV, then 3.16 TeV < 
Mz" < 6.78 TeV. For a different M s the range of Mz» is altered because of changes in g[, 9, and 0; e.g. 
for M s = 10 19 GeV, the range becomes 2.82 TeV < M z „ < 6.05 TeV. 
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FIG. 3: An exhibition of the SM ++ vacuum stability patterns in the rrih" vs a plane, for fixed 
masses of Z"; M z » = 3.5 TeV (left pannels), M z » = 4.5 TeV (middle pannels), and M z » = 6 TeV 
(right pannels). The analysis is based on a scan of 10 8 trial random points with M s = 10 14 GeV 
(upper row), M s = 10 16 GeV (middle row), and M s = 10 19 GeV (lower row). The points yielding 
a stable vacuum solution up to M s are blue-printed, those leading to unstable vacuum solutions 
are red-printed, and points giving runaway solutions (i.e., those in which the Higgs doublet self- 
coupling blows up) are purple-printed. 
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FIG. 5: Typical behavior of the running couplings A,(t), for the average value of the initial condi- 
tion, (Aj° ) = 0.25 in the integration of (25). 
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FIG. 6: The allowed SM ++ parameter space in the mh" vs a plane under the vacuum stability 
constraint of Eq. (5), for the case M z » = 4.5 TeV, with M s = 10 14 GeV (blue), M s = 10 16 GeV 
(green), and M s = 10 19 GeV (red). The perturbativity upper bound is defined by Aj < 27T. 



solutions are red-printed, and points giving runaway solutions are purple-printed. A 
stable vacuum solution is one in which the positivity condition (5) is fulfiled all the 
way to A = M s . An unstable solution is one in which the stability conditions of the 
vacuum (Ai > 0, A 2 > 0, AiA 2 > A§/4) are violated. A runaway solution is one in 
which the RG equations drive the Higgs doublet self-coupling non-perturbative. The 
perturbativity upper bound (sometimes referred to as 'triviality' bound) is given by 
Ai < 27r at any point in the RG evolution [20] . In the first row of Fig. 3 we display 
sections of constant Mz» for the scan shown in Fig. 2. In the second and third rows 
of Fig. 3 we show theree different sections of the myi — a plane, corresponding to the 
same values of Mz», but for scans with different string scales: M s = 10 16 GeV (second 
row) and M s = 10 19 GeV (third row). The vacuum stability condition is driven by the 

behavior of Ai, and actually is largely dominated by the initial condition A^. Indeed, 
if the extra gauge boson Z" gets its mass through a non-Higgs mechanism and the 
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FIG. 7: Variation of SM ++ vacuum stability regions with M Z ". We have taken M s = 10 16 GeV, 
M z » = 3.5 TeV (red), M { Z") = 4.5 TeV (green), and M(Z") = 6.0 TeV (blue). The perturbativity 
upper bound is defined by \ < 2tt. 
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FIG. 8: Variation of SM ++ vacuum stability regions with m^. We have taken a = 0.06, Mz" = 
4.5 TeV, M s = 10 14 GeV (blue), M s = 10 16 GeV (green), M s = 10 19 GeV (red). The perturbativity 
upper bound is defined by Aj < 2tt. 

scalar potential (3) is that of SM (i.e. v" = \ 2 = A3 = 0), the RG evolution collapses 
to that of SM and there are no stable solutions. 3 



3 Of course, even if v' 1 



A2 — An 



0, with an extra gauge boson the RG evolution of Ai is not exactly 
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FIG. 9: Variation of SM ++ vacuum stability regions with g'^M s ). The stable regions correspond 
to g[(M s ) = 1.0 (cyan), 0.85 (blue), 0.6 (green), 0.485 (red), 0.4845 (purple). We have taken 
M s = 10 14 GeV, M Z " = 4.5 TeV, m h = 125 GeV. The perturbativity upper bound is defined by 
Xi < 2tt. 

• To determine the range of initial conditions A^ yielding stable vacuum solutions we 
fit the boundaries of the blue band in the scatter plot, with M s = 10 14 GeV. The 
results shown in the left panel of Fig. 4 lead to 0.19 < A^ < 0.36. The lower limit 
of which defines the boundary between stable and unstable solutions, is close to 
the value required for vacuum stability of the SM potential, as shown in (1). Namely, 
substituting rrih = 130 GeV and a = in (13) we obtain Aj° = 0.14. The similarities 
between the minimum value of rrih that allows absolute stability up to the Planck scale 
within SM and the minimum value of rrih in the decoupling limit of (13) reinforces 
our previous statement concerning the strong dependence of the RG evolution with 
the initial condition \[ . The allowed range of initial conditions for stable vacuum 
solutions depends on the value of the string scale; e.g. for M s = 10 19 GeV we obtain 
0.22 < A[ 0) < 0.30. We have also determined the average value of the initial condition 
Aj through a fit to the blue points in the scattered plot. The results shown in the right 
panel of Fig. 4 lead to (A 1 °' ) ) = 0.25, independently of the value of M s . The typical 
behavior of Ai, A2, and A1A2 — A§/4, for the average value of the intial condition (A^), 
is display in Fig. 5. 

• To determine the sensitivity of the RG evolution with respect to the triviality bound we 
duplicate the analysis, but setting the perturbativity upper bound in all the couplings: 
Xi < 2tt. The contours display in Fig. 6 (for Mz" = 4.5 TeV) show there is almost no 
visible departure from the contours in Fig. 3. 



that of SM, see (25). 
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• In Fig. 7 we display the sensitivity of the RG evolution with M z ». For large values 
of a there is no variation in the contour regions. For a < 0.03 there are some small 
variances. This small differences show the effect of the initial conditions A^ and Ag ^ 
(both depending directly on M z ») on the evolution of the system. 

• We have verified that there is no significant variation of the SM ++ vacuum stability 
regions within the uncertainty. An example for a = 0.06 and M Z " = 4.5 TeV is 
given in Fig. 8. 

• Finally, in Fig. 9 we display the variation of our results with g[(M s ). It is clearly seen 
that for g[(M s ) > 0.4845 the dependence on g[ seems to be fairly weak. However, 
if g[(M s ) < 0.4845, the g' 4 coupling becomes non-perturbative and the stable region 
of parameters shrink dramatically. The stability of SM ++ vacuum is then nearly 
independent of the Z" branching fractions [24]. 

• The low energy effective theory discussed in this paper requires a high level of fine 
tuning, which is satisfyingly resolved by applying the anthropic landscape of string the- 
ory [35-37]. Alternatively, the fine tuning can be circumvented with a more complete 
broken SUSY framework. Since in pure SUSY the vacuum is automatically stable, 
the stability analysis perforce involves the soft SUSY-breaking sector. Hence rather 
than simply searching for the Higgs self-coupling going negative in the ultraviolet, the 
stability analysis would involve finding the local and global minima of the effective 
potential in the multi-dimensional space of the soft-breaking sector [38]. However, the 
Higgs mass range favored by recent LHC data may be indicative of high-scale SUSY 
breaking [39]; perhaps near the high energy cutoff of the field theory, beyond which a 
string description becomes a necessity [40]. 

• In summary, we have shown that SM ++ is a viable low energy effective theory, with 
a well-defined range of free parameters. It is of immediate interest to calculate the h 
couplings (after diagonalization) to W + W~ , bb, tt and determine whether there is a 
region of the parameter space that would allow an enhancement in h — > 77, as hinted 
in LHC data [41]. 
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